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Abstract. We discuss numerical and theoretical results for models of magnetization 
switching in nanoparticles and ultrathin films. The models and computational methods 
include kinetic Ising and classical Heisenberg models of highly anisotropic magnets 
which are simulated by dynamic Monte Carlo methods, and micromagnetics models of 
continuum-spin systems that are studied by finite-temperature Langevin simulations. 
The theoretical analysis builds on the fact that a magnetic particle or film that is 
magnetized in a direction antiparallel to the applied field is in a metastable state. 
Nucleation theory is therefore used to analyze magnetization reversal as the decay of 
this metastable phase to equilibrium. We present numerical results on magnetization 
reversal in models of nanoparticles and films, and on hysteresis in magnets driven by 
oscillating external fields. 

1 Introduction 

In recent years, the interest in nanostructured magnetic materials has soared for 
a variety of reasons. For one, it is only quite recently that it has become possible 
to synthesize and measure nanometer-sized magnetic particles in small, ordered 
arrays, often by techniques that involve modern, atomic-resolution microscopies, 
such as scanning-tunneling microscopy (STM), atomic force microscopy (AFM), 
or magnetic force microscopy (MFM) . AFM and MFM pictures of 

an array of nanometer-sized iron pillars, fabricated by STM-assisted chemical 
vapor deposition are shown in Fig. |l|. The techniques are currently becom- 
ing precise enough to even allow investigation of individual nanoparticles. At 
the same time, computers have had a profound influence in two different ways. 
The need for ever higher data-recording densities has driven the size of parti- 
cles used in recording media down into the nanometer range 0,^,^, while the 
rapidly increasing power of computers has made it feasible to perform simula- 
tions of the dynamic properties of realistic model systems of sizes comparable 
to experimental ones |l0|[ll| , p^ . 
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Fig. 1. Array of nanoscopic iron pillars of dimensions approximately 40 x 40 x 200 nm^, 
grown by STM-assisted chemical vapor deposition, (a) AFM image of the array of pil- 
lars, as grown on top of a ^m-size Hall-effect magnetometer, (b) MFM image of the 
array after thermal randomization in near-zero applied field. The magnetic field from 
each pillar is imaged and seen to point along the major axis of the pillars, either up 
(white) or down (black), (c) MFM image of the array in an applied field of 200 G. Al- 
most all the magnets are aligned with the field. Image data courtesy of D.D. Awschalom 
and J. Shi. (a) and (b) a,fter |^ A pdf version of this paper, which includes this fig ure, 
can be downloaded from tittp://www/csit.fsu.edu/'~ rikvold/abstracts/turkey.pdl 



Modern magnetic recording technologies involve particles that are near the 
superparamagnetic limit. In this limit, the energy barrier separating the two en- 
ergetically degenerate magnetic orientations is small enough that thermal fluc- 
tuations frequently lead to spontaneous switching of the orientation. As a result, 
the magnetic coercivity decreases with decreasing particle size for particles be- 
low the superparamagnetic size limit. The limit appears as a maximum in a 
curve showing switching field or coercivity versus particle size, such as in Fig. |[ 
Since the random magnetization reversals in particles below the superparamag- 
netic limit degrade recorded information, the engineering challenge has been to 
keep the energy barrier in the individual particles high enough to make sponta- 
neous switching infrequent while keeping the material magnetically soft enough 
to facilitate recording. As the volumes of the magnetic particles have shrunk 
to reach recording densities on the order of 100 Gb/in^ or more mate- 
rials with higher cocrcivitics due to strong crystalline anisotropics have been 
employed iQ. In order to enhance engineering practices, it is essential to extend 
the physical understanding of the superparamagnetic limit past the theories of 
uniformly mangetized particles to include magnetization reversal dynamics that 
proceed through localized regions of reversed magnetization that subsequently 
spread throughout the magnetic element. 

The remainder of this article is organized as follows. In Sec. 2 we summarize 
some aspects of the theory of magnetization switching in anisotropic magnets, 
including effects of anisotropy (Sec. 2.1), nucleation theory (Sec. 2.2), and model 
systems (Sec. 2.3). In Sec. 3 we give some new results of finite-temperature micro- 
magnetic simulations of magnetic nanoparticles. In Sec. 4 we discuss hysteresis in 
nanoparticles and ultrathin films, in particular the frequency dependence of hys- 
teresis loops (Sec. 4.1) and a dynamic order-disorder phase transition (Sec. 4.2). 
A brief summary and conclusions are given in Sec. 5. 

2 Theory of Magnetization Switching in Anisotropic 
Magnets 




Fig. 2. Effective switching field (analogous to coercivity) versus particle size for single- 
domain ferromagnetic barium ferrite particles (o with error bars, right vertical and top 
horizontal axes, experimental results) and for two-dimensional L x L Ising systems 
(data points connected by solid lines, left vertical and bottom horizontal axes). The 
barium ferrite data were digitized from Fig. 5 of [^^. The Ising data are Monte Carlo 
simulations from for waiting times r=100 Monte Carlo Steps per Site (MCSS) (o) 
and 1000 MCSS (x) at r=0.8Tc, where Tc is the exact Ising critical temperature. After 



2.1 Effects of Magnetic Anisotropy 

The most common description of magnetization switching is the mean-field, 
uniform- rotation theory of Neel ||l^ and Brown ||l^,|2^. One assumes uniform 
rotation of all localized moments in the particle to avoid an energy barrier due 
to exchange interactions of strength J. The remaining barrier, Z\, is caused by 
magnetic anisotropy - a combination of crystal-field and magnetostatic effects. 
The equilibrium thickness of a wall separating oppositely magnetized domains is 
£w -y/ J/ A. For particles smaller than with small anisotropy, the uniform- 
rotation picture is reasonable. If the anisotropy is largely magnetostatic, the 
competition between exchange interactions and the demagnetizing field favors 
domains of opposite magnetization in particles larger than The domains 
control switching through the field-driven motion of preexisting domain walls 
[ pTp^ , ^ . However, if the anisotropy is largely crystalline, there exists a range 
of single-domain particle sizes that are larger than but smaller than the size 
at which the particle becomes multidomain (often the case in ultrathin films 

IE)- 
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In anisotropic nanomagnets the state of uniform magnetization opposite to 
the apphed field constitutes a metastable phase. This nonequilibrium phase de- 
cays by thermaUy assisted nucleation and subsequent growth of locahzed regions, 
inside which the magnetization is parallel with the field |]T6j . These growing re- 
gions are referred to as droplets to distinguish them from equilibrium domains. 
This mechanism yields results very similar to recent experiments on single- 



domain nanoscale ferromagnets |25 



2.2 Application of Nucleation Theory to Magnetization Reversal 

Here we present a short summary of homogeneous nucleation theory as it applies 
to uniaxial magnets. This theory covers situations in which the switching events 
are nucleated by thermal fluctuations, without the influence of defects. Further 



details are available in |lX|,|lJ,|17|j26||23,|28|,g^ 

The central problems in nucleation theory are to identify the fluctuations 
that lead to the decay of the metastable phase and to obtain their free-energy 
cost relative to the metastable phase. For anisotropic systems dominated by 
short-range interactions, these fluctuations are compact droplets of radius R. 
The free energy of the droplet has two competing terms: a positive surface term 
cx a{T)R'^~^ and a negative bulk term cx \H\R'^ where d is the spatial dimension, 
a{T) is the surface tension of the droplet wall, and H is the applied magnetic field 
along the easy axis. Their competition yields a critical droplet radius, Rc{H, T) cx 
a{T)/\H\. Droplets with R < Rc most likely decay, whereas droplets with R > 
Rc most likely grow to complete the switching process. The free-energy cost of 
the critical droplet (i? = Rc) is AF{H, T) oc a{TY /\H\'^-'^ . Nucleation of critical 
droplets at nonzero temperature T is a stochastic process with nucleation rate 
per unit volume given by an Arrhenius relation: 

I{H, T) cc \H\^ exp [-PAFiH, T)] \H\^ exp [-pSiT)/\H\^-'] , (1) 

where /3 = l/fceT (fee is Boltzmann's constant), S{T) is the i7-independent 
part of AF, and the prefactor exponent K is known for many models from field- 
theoretical arguments ||l7| , ^ , ^Tl|3^ . The particles are of finite size L, and the 
dominant reversal mechanism depends on H, T, and L. 

In the weakest applied fields, the particles are in the "Coexistence" (CE) 
regime, with the average metastable lifetime tce{H, T, L) ~ exp \2l3a{T)L'^~^^ . 
(This result is nearly independent of the boundary conditions |27| .) The regime 
corresponds to R^ > L, and the associated L-dependent crossover field is called 
the Thermodynamic Spinodal (ThSp) [p^ , p8[|29| . Estimating its value by assum- 
ing Rc{H, T, L) K, L, one finds HThSuiT, L) ^ L^^. The L-dependence of -ffThSp 
is given by the dotted curve in Fig. g. 

For \H\ > H-ThSp (but not too large), the lifetime is determined by the inverse 
of the total nucleation rate, 

TsB{H,T,L)^{L''l{H,T)y' ^L-''\H\^exp[(3S{T)\Hf-^] , (2) 
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Fig. 3. Lifetimes for two-dimensional I?" Ising systems with i = 20 and 160 at 
T — 0.57Tc (a), and three-dimensional Ising systems with L = 10 and 40 at 
T = OAATc (b). The data points are direct Monte Carlo simulation results, while 
the lines are extrapolations with the Projected Dynamics (PD) accelerated dynam- 
ics algorithm p3| , p^j35[ , based on the smallest system at the weakest field. The sharp 
changes in slope correspond to the DSp, with the deterministic regime to the left and 
the stochastic regime to the right. The ratio of the slopes of the curve in the single- 
droplet and multidroplet regimes is (d+l), in agreement with and Q. After ^] 



which is inversely proportional to the particle volume, L'^ (see Fig. The 
subscript SD stands for Single Droplet and indicates that in this regime the 
switching is normally completed by the first droplet to reach R^. In both of 
the stochastic reversal regimes(CE and SD) the probability that switching has 
not taken place within a time t after the field reversal, Pnot{t), takes the form 
PnotW =exp(-t/r). 

A second crossover, called the Dynamic Spinodal (DSp) [ pTp^ , ^ , is a con- 
sequence of the finite velocity, v w of the surface of a growing supercrit- 
ical droplet A reasonable criterion to locate the DSp is that the average 
time between nucleation events, which is tsd, should equal the time it takes a 
droplet to grow to a size comparable to L. This yields the asymptotic relation 
Husp{T,L) ~ nn(L)]i/('^-i). The L-dependence of ^dSp is given by the dashed 
curve in Fig. g. For \H\ > 7?DSp, the metastable phase decays through many 
droplets which nucleate and grow independently in different parts of the system. 
In this Multidroplet (MD) regime [oj2a,B9], the classical Kolmogorov-Johnson- 



Mehl-Avrami (KJMA) theory of metastable decay in large systems [37, 38| , |39[ 
gives the lifetime 

tmb{H,T) ^ [l{H,T){iy\H\)^/{d+l)ln2y'^''''^'^ , (3) 

independent of L (see Fig. |^). In the MD regime Pnot(0 ~ erfc {{t — t)/A) p^ , 
where the width A of the switching-time distribution depends on iJ, T, and L. 
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Fig. 4. Metastable phase dia- 
gram for tiie two-dimensional 
Ising model at temperatures up 
to Tc. The dashed curve rep- 
resents HuFSpiT) ■ Data points 
connected by solid line seg- 
ments represent HoSp{T) for 
several values of L between 
20 and 10^ The filled data 
points are the results of Monte 
Carlo simulations, while the 
empty data points represent 
a low-temperature approxima- 
tion [M for large systems 



For very strong fields nucleation theory becomes irrelevant to the switching 
behavior. A reasonable way to estimate the crossover to this Strong-Field (SF) 
regime is to require that the critical radius should be on the order of the lattice 
constant a. Specifically, requiring Rc = a/2, we get the crossover field called the 
mean-field spinodal (MFSp) HuFSpiT) « 2cr(T)/meq(r), where moq(T) is the 
zero-field equilibrium magnetization. The "metastable phase diagram" in Fig. || 
shows Hmfsp{T) , as well as HDSp{T) for two-dimensional Ising systems of widely 
varying sizes. Note the logarithmically slow convergence to zero of -ffosp with 
increasing system size. As a result, even macroscopic metastable systems may 
be "small" in the sense that they decay via the single-droplet mechanism. 

The switching field, 77sw(iw, T, L), is the field required to observe a specified 
average waiting time, t^. It is found by solving r in the relevant region (CE, SD, 
or MD) for H with T=t^. The resulting L dependence of i?sw is a steep increase 
with L in the CE regime, peaking near the ThSp, followed by a decrease in 
the SD regime towards a plateau in the MD regime [[l6|j26| . This behavior is 
illustrated in Fig. |^. Note that the maximum in Hgv, (related to the maximum 
coercivity) occurs even in the absence of dipole-dipole interactions. For other 
boundary conditions and in systems with dipole interactions the i?sw versus L 
curve can even have more that one maximum |p7| . Maximizing the coercive field 
is important in magnetic recording applications. 

2.3 Statistical-mechanical Model Systems 

Simplified statistical-mechanical models are often amenable to analytic solu- 
tions, with no free fitting parameters, that agree well with the numerical results. 
Despite their lack of realism, they are therefore important as testing grounds 
for theoretical descriptions of different switching mechanisms. In more realistic 
models, that correspond to larger numbers of actual materials, analytic results 
are difficult to come by, but the physical insights gained from the simpler mod- 
els can be readily applied. Here we introduce three such models in order of 
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increasing complexity: the kinetic Ising model, the classical Heisenberg model, 
and finite-temperature Langevin micromagnetic models. 



The Ising Model The simplest microscopic model of a ferromagnet is the 
nearest-neighbor Ising model, in which discrete spins, Si = ±1, are placed on 
the sites (labeled i) of a two- or three-dimensional lattice. The spins interact 
with their neighbors with a strength J, so that the model is described by the 
Hamiltonian 

H=-JJ2''^^J-HJ2s^. (4) 

The model can easily be generalized to longer-range interactions, different lattice 
geometries, etc. Despite its apparent simplicity, it has many of the attributes of 
more complicated systems, while many of its properties are exactly known. It is 
therefore a very commonly studied model. 

The Ising Hamiltonian, (^), is not a true quantum- mechanical Hamiltonian, 
and the Ising model therefore docs not have an intrinsic dynamic. To simulate 
thermal fluctuations one uses Monte Carlo simulation of a local stochastic dy- 
namic which does not conserve the order parameter. An often-used example is 
the Metropolis dynamic with the spin-flip probability 

WmWAE) = min [1, cxp {-/SAE)] , (5) 

where AE is the energy change that would ensue if the flip were to occur. 
Another popular choice is the Glauber dynamic defined by 

The basic time scale of the Monte Carlo simulation is not known from first 
principles, but it is expected to be on the order of a typical inverse phonon 
frequency, 10~^-10~^^ s. In dynamics such as these, where each potential flip 
is accepted or rejected randomly, flips can become very rare when rejection 
rates are high. To perform simulations on the very long time scales necessary to 
observe metastable decay, one needs to use rejection-free Monte Carlo algorithms 



5|,[46[ and other advanced algorithms |^ 34,35||. 
Analogous stochastic time evolutions can also be imposed on models whose 
spins have continuous degrees of freedom. Here we briefly discuss one such model, 
the anisotropic Heisenberg model. 



The Anisotropic Heisenberg Model Like the Ising model, the Heisenberg 
model consists of spins located at discrete points on a lattice. However, unlike 
the spins in the Ising model, which equal ±1, Heisenberg spins are n-dimensional 
vectors of unit length. When n = 2, this model is usually referred to as the XY or 
plane-rotor model. The Hamiltonian for the nearest-neighbor Heisenberg model 
with only interaction anisotropy is 

H = - ^ ^ [J„s„,iS„j] - if ■ ^ Si , (7) 
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where Sn^i is the n-th component of the i-th spm vector s^, J„ are coupUng 
constants, and H is the external magnetic field. For the example of this model 
presented here, n = 3, Jx — Jy — ^, Jz = '2, H ~ HzZ, and the lattice is a 
two-dimensional L x L square lattice. 

There are many stochastic dynamics for the Heisenberg model which yield 
identical equilibrium results, but have different relaxation dynamics. The dy- 
namic assumed here consists of selecting a spin vector at random, then choosing 
a new orientation for that spin, uniformly distributed over the unit sphere, and 
then accepting or rejecting the new configuration based on (^. Other dynam- 
ics exist which make only small changes to spin orientations, but these are not 
discussed here. 
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Fig. 5. Switching behavior for an anisotropic n = 3 Heisenberg model on a two- 
dimensional L X L square lattice. The parameters are Jx = Jy = 1, Jz = 2, and T = 1, 
which is below the critical temperature for this model. The circles are for a system with 
L = 16, and the single-droplet (SD), multidroplet (MD), and strong-field (SF) regimes 
are labeled. The insets show typical system configurations during the switching process 
for systems with L — 64. Note that the dynamic spinodal (DSp) depends on L, and 
thus the SD and MD regimes appear at different fields for L = 16 and L = 64. The 
grayscale for the insets shows with lighter shades indicating metastable spins and 
darker shades indicating more stable spins 
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The simulation begins with all Si = —z and Hz > 0. This mctastable phase 
then decays in a manner consistent with homogeneous nucleation and growth. 
However, unlike the Ising model, the continuous degrees of freedom add addi- 
tional complications, such as effective long-range interactions between droplets. 
We have not yet attempted to quantify these differences. Figure || shows lifetimes 
and configuration snapshots in the single-droplet, multidroplet, and strong- field 
regimes for an anisotropic Heisenberg model at a temperature below criticality. 
The field dependence of the lifetime is seen to be very similar to that of the 
two-dimensional Ising model, shown in Fig. ^(a). 



M. 




Fig. 6. Magnetization along 
the pillar long axis, Mz, at 
three times during the switch- 
ing process. Light shades rep- 
resent the metastable orienta- 
tion and dark shades the equi- 
librium orientation. The micro- 
magnetic simulation shown in 
this figure models an individ- 
ual pillar using a 7 x 7 x 49 
lattice at T = 20 K. The pillar 
is shown in a one-quarter cut- 
away view 



Finite-temperature Langevin Micromagnetics More realistic representa- 
tions of nanoscale magnetic systems can be obtained by micromagnetic modeling. 
In this method the "spins" are coarse-grained magnetization vectors M{ri); each 
represents the magnetization within a cell centered at position r^. In this low- 
temperature model ^7| , the vectors have a fixed magnitude Mg corresponding 
to the bulk saturation magnetization density. The time evolution of each spin 
is governed by the damped processional motion given by the Landau-Lifshitz- 
Gilbert (LLG) equation ||,|| 

dM{r,) 



dt 



^^M(rO X (H{r,) - ^M(rO x H{r., 



(8) 
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where the electron gyromagnetic ratio is 70 = 1.76 x 10^ Hz/Oe and a is 
a phenomenological damping parameter. The local field at the i-th spin, H{ri), 
is generally different at each location. This field is a linear superposition of 
fields, one for each type of interaction in the system. Typical examples include 
fields from external sources, exchange, crystalline anisotropy, and dipole-dipole 
interactions. Thermal fluctuations may also contribute a term: a stochastic field 
i?n(^i) that is assumed to fluctuate independently for each spin |2^. The fluc- 
tuations are assumed Gaussian, each with zero first moment and with the second 
moments given by the fiuctuation-dissipation relation po| 



where ffn^ indicates one of the Cartesian components of H^. Here V — (Ar)^ 
is the discretization volume of the cell, is the Kronecker delta representing 
the orthogonality of the Cartesian coordinates, and S{t — t') is the Dirac delta 
function. 

While this stochastic term necessitates careful treatment of the numerical 
integration in time of this stochastic differential (Langevin) equation, the most 
computationally intensive part of the calculation involves the dipole-dipole term. 
For systems with more than a few hundred model spins, it is necessary to use 
a sophisticated algorithm such as the fast multipole method (FMM) ||5^,|l[ . An 
extensive discussion of the issues involved in finite-temperature simulations of 
micromagnetics is presented in ||l^ . The growth of a droplet during the switching 
of an iron nanopillar at T=20K is shown in Fig. ||. Nucleation is observed to 
occur at the ends of the pillars p^ . 

3 Finite-temperature Micromagnetics Results for 
Nanoparticles 

Micromagnetic simulations have been applied to the study of iron pillars mod- 
eled after those shown in Fig. |l|. Unless otherwise noted, the model iron pillars 
discussed here are 5.2 nm x 5.2 nm x 88.4 nm. The cross-sectional dimensions 
are small enough, about two exchange lengths, that the only significant inhomo- 
geneities in the magnetization are those in the z direction, i.e. along the long 



axis of the pillar [p2y53| . In light of this, the pillars have been modeled as a linear 
system of magnetic cubes with side 5.2 nm. This model, discussed previously in 
[ |l2|j54[] , includes thermal fluctuations, exchange, and dipole-dipole interactions. 

The results for Pnot{t), with r=100K are shown in Fig. for applied fields 
of iJ=1000Oc and 800 Oe. Here switching is defined to occur when the z- 
component of the total magnetization, M^, passes through zero. The form of 
Pnotit) is not exponential, which can be explained by the fact that nucleation of 
the reversed droplets is easier at the ends of the pillars than in the middle. As- 
suming that the nucleation rate for droplets at the end of the pillars is constant, 
/, and that the earliest time switching can occur because of the finite velocity 
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of droplet growth is to, it can be shown that the probability of not switching is 



The parameters are fitted by matching the first and second moments of the 
simulation results to those of the theoretical forms. As long as the applied field 
is relatively weak, the agreement between the theory and the simulations is quite 
good. Switching at times to < t < 2to is possible only when nucleation occurs at 
both ends. 

Since the nucleation occurs at the ends, the dependence of the switching on 
the size of the system is different from that seen in isotropic models. Results for 
the parameters to (squares) and l/I (circles) at T=100K and _ff=1000Oe are 
shown in Fig. ^ for pillars of different lengths, i.e. composed of different numbers 
of cubes 5.2 nm on a side. The nucleation rate is nearly constant, indicating that 
the size of the energy barrier does not depend on the pillar length. The growth 
time, indicated by to, however, increases as the droplets have to grow farther to 
switch the magnetization. The nearly linear increase with pillar length indicates 
that the interface velocity is not significantly affected by the demagnetizing field 
associated with the high aspect ratio of the pillars. 

Finally, changes in the switching mode as the field is changed are shown 
in Fig. p. Here the mean switching time, (tsw), and standard deviation, at, 
are shown versus applied field for the 88.4 nm long pillars at r=100K. The 
mean, (tsw), and standard deviation, at, of the switching time tsw versus inverse 
applied field for pillars of the same type as considered in Fig. ^ At weak fields 
the mean and standard deviation are nearly equal, where the exponential tail at 
t > 2to dominates (|To|). As the applied field is increased, the barrier to nucleation 




t < to 

to<t< 2to 
2to < t 



(10) 








Fig. 7. Probability of not 
switching before time t, 
Pnotit), for micromagnetic 
simulations at r=100K 
with applied fields of 
1000 Oe and 800 Oe as 
labeled. The solid curves 
are simulation data for 624 
and 252 switches, respec- 
tively, for pillars modeled 
as a one-dimensional chain 
of spins. The dashed curves 
are fits to the theoretical 
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decreases, and the exponential behavior becomes less dominant. Eventually, ( fi^ 
breaks down as the multidroplet reversal mechanism becomes important. The 
multidroplet nature of the reversal has been verified by direct observation of the 
switching, which shows droplets nucleating away from the ends at i?=1000Oe. 



4 Hysteresis 

Hysteresis is common in many nonlinear systems driven by an oscillating external 
force, including nanostructured magnets in an oscillating field. It occurs when 
the dynamics of the system is too sluggish to keep pace with the force. The term 
was coined by Ewing in the context of magnetoelasticity from the Greek 
word husterein {varepeuj) which means "to be behind." 



4.1 Hysteresis- loop Areas 

Among the earliest aspects of hysteresis to receive sustained interest is the 
hysteresis-loop area. In the magnetic context of this article, the hysteresis loop is 
a plot of magnetization versus applied field, and its area is given by the integral 
A — — § m{H)<lH . A typical hysteresis loop for a small system with thermal 
noise is shown in Fig. |lO|. The particular importance of the loop area is that it 
corresponds to the energy dissipation per period of the applied field. It is thus 
relevant to the performance of most electrical and electronic equipment. Recent 
experiments on ultrathin Fe and Co films with Ising-like anisotropy have con- 
sidered the frequency dependence of the hysteresis-loop areas [011,11,13. The 
results of these studies were interpreted in terms of power laws, but with expo- 
nents that vary widely between experiments. The experimental situation thus 
may appear somewhat unclear. 
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Fig. 8. Inverse nucleation 
rate, /~^, and earliest 
switching time, to, from 
(p^), for different pillar 
lengths. The rate of nucle- 
ation at the ends of the pil- 
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on the pillar length, and 
therefore likewise the en- 
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Fig. 10. Hysteresis 
loop for one model Fe 
pillar of the same type 
considered in Fig. ^ 
at T = 100 K for an 
oscillating field with a 
period of 75 ns 



A resolution is provided by the nucleation-and-growth picture of magne- 
tization switching presented here. We assume that the system is driven by 
a sinusoidally oscillating field, H{t) — Hosin^Lut). Since the nucleation rate, 
which has the dimension of a frequency per unit volume, is proportional to 
exp [— /3S'(T)/|if|'^~^] by (p, one would expect that the field at which the mag- 
netization changes sign should depend on the frequency as — (Ino;)"^''*-''"^-'. The 
loop area is approximately proportional to the switching field multiplied by the 
saturation magnetization (see Fig. Thus, one would expect the loop area 
to show this logarithmic frequency dependence in the asymptotic low-frequency 
limit. Analytic calculations have confirmed this asymptotic result. However, for 
higher, but still low, frequencies they show a very slow crossover to the asymp- 
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iogio(i/R) 

Fig. 11. The hysteresis-loop area, A, versus the dimensionless frequency, 1/R. (Top) 
Two-dimensional Ising model of an ultrathin film. Data points: Monte Carlo simulations 
for L = 64, r=0.8Tc, and Ho=0.3J. For these parameters the magnetization switch- 
ing occurs via the multidroplet mechanism, except for the lowest frequencies. Solid 
curve: numerical integration for sinusoidally varying field. Dotted curve: numerical in- 
tegration for linearly varying field. Dot-dashed curve: numerical integration assuming 
magnetization reversal via the single-droplet mechanism. Dashed curve: low-frequency 
asymptotic solution. Power-law fits would yield very different effective exponents for 
fits centered at different frequencies. After jHoj. (Bottom) Micromagnetic model of 
iron pillars with length 88.4 nm and square cross-section 5.2 nm x 5.2 nm for r=100K 
(circles) and r=20K (squares) with _ffo=2000 Oe. The lines are guides to the eye. The 
inset shows the correlation between the magnetization and field versus frequency. The 
lowest-frequency zero-crossing, indicating a resonance condition, occurs at roughly the 
same frequency as the maximum in A 



totic behavior, which is confirmed by Monte Carlo simulations. Such a slow 
crossover could easily be mistaken for a power law, even when observed over 
several decades in frequency |5fl,pl^,p2|. This result was shown to hold, both 
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when the magnetization reversal occurs via the single-droplet mechanism p2[ 
and the multidroplet mechanism p6[ , even though the details are different. The 
behavior is illustrated in the top part of Fig. |l^. Recently we have also found 
analogous behavior in micromagnetics simulations of nanometer-sized iron pil- 
lars, see the bottom part of Fig. |ll|. In these figures the frequency is given in 
terms of the dimensionless frequency 1/i? = ujt{Hq, T)/2tt. 



4.2 Dynamic Phase Transition 

Different phenomena occur in hysteretic systems as the driving frequency is in- 
creased. Eventually the field will vary too quickly for the system to have time to 
switch during a single period. Thus, low frequencies lead to symmetric hysteresis 
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Fig. 12. The dynamic phase transition in a two-dimensional Ising system at T = O.STc, 
driven by a square- wave oscillating field of amplitude Hq = 0.3 J. After (a) The 
dynamic order parameter, (IQI), versus the dimensionless period 7? for several system 
sizes, (b) The order-parameter fluctuation strength, X^, versus R for several system 
sizes 



loops, such as in Fig. while high frequencies produce loops in which the mag- 
netization oscillates about one or the other of its zero- field equilibrium values. 
For small systems or weak field amplitudes, such that the magnetization reversal 
occurs via the single-droplet mechanism, this results in stochastic resonance ||6^ . 

For large systems or stronger fields, such that the magnetization switching 
occurs via the multidroplet mechanism, the transition from symmetric to asym- 
metric hysteresis loops becomes a genuine critical phenomenon at a sharply 
defined critical frequency. The transition is essentially due to a competition be- 
tween two time scales: the metastable lifetime, t{Hq, T), and the frequency of the 
applied field, lo/2ti. As a result, the critical value of the reduced frequency is 
on the order of unity. This nonequilibrium phase transition was first observed in 
numerical solutions of mean-field equations of motion for ferromagnets in oscil- 
lating fields 1^,0 • Subsequently it has been observed in numerous Monte Carlo 
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simulations of kinetic Ising systems J5|f|,|6|j67[||,|9|,^o|j75,|7|,||,0|7| and in 
further mean- field studies ||,^0|j75,|73|,|76f. It may also have been experimentally 
observed in ultrathin films of Co on Cu(lOO) [|§7^. 

In this far-from-equilibrium phase transition the role of order parameter is 
played by the period-averaged magnetization, Q ~ [u j^ii) § m(t)dt. This quan- 



tity is shown in Fig. 12(a) versus the dimensionless period R for several system 
sizes. The order-parameter fluctuation strength, = [{Q^) — {\Q\)^] j which 
corresponds to the susceptibility in an equilibrium system, is shown for several 
values of L in Fig. |l^(b). Both the order parameter and its fluctuations de- 
pend on L in a way very similar to data from simulations of equilibrium phase 
transitions. And, indeed, formal finite-size scaling analysis of the Monte Carlo 
data ||5^ , |65| , [75| , as well as analytical arguments [fz 8(1791 , have shown that this 
far-from-equilibrium phase transition belongs to the same universality class as 
the equilibrium phase transition in the Ising model in zero field. This is a quite 
remarkable result, as it extends the scope of an equilibrium universality class to 
a far-from equilibrium system. 



5 Summary 

In this article we have presented numerical and theoretical results on magne- 
tization reversal and hysteresis in models of magnetic nanoparticles and ultra- 
thin films. Models that were explicitly considered are kinetic Ising and clas- 
sical Heisenberg models, which were studied by dynamic Monte Carlo sim- 
ulations, and continuum-spin micromagnetics models, which were studied by 
finite-temperature Langevin-equation methods. The simulation results were in- 
terpreted within the context of nucleation theory, and it was shown how the 
reversal modes change from single-droplet to multidroplet upon increasing the 
strength of the applied field or the size of the system. Computer simulations of 
model systems such as those presented here enable one to study in detail the 
statistical properties of the reversal processes, as well as the time dependent 
internal magnetization structure. Such simulation results have now attained suf- 
ficient quality that they can fruitfully be compared with present and future 
experiments. 
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